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Abstract 

A perpendicular ion drift is investigated as a mechanism for the 
generation of magnetic field structures in a highly collisional dusty 
plasma. The basic dissipation mechanism is assumed to be the dust- 
neutrals momentum exchange, so that plasmas with a small ionization 
fraction are natural candidates for experiments. The model reduces 
to a nonlinear partial differential equation for the vector potential. 
The conditions for linear instability are determined. The stationary 
states are periodic arrangements for the magnetic field, described by a 
Lienard equation. The fully depleted (ion-dust) case is also considered 
in detail. Applications to planetary rings, comets and dust acoustic 
instability experiments are analyzed. A necessary condition for the 
validity of the model is a sufficiently slow time-scale of the generated 
magnetic fields. 

1 Introduction 

Collisional processes have an important role in dusty plasmas in laboratory 
or in the environment, as for the dust-acoustic and dust-ion-acoustic wave 
propagation and instabilities [HE]- In particular [3], collisions between the 
charged dust particles and the neutrals are relevant, in view of the usually 
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small degree of ionization in applications. Ion-neutrals collisions are also 
fundamental in experimental observation of dust-acoustic waves [I], while 
dust-neutrals collisions are related to the generation of dust-acoustic waves 
driven by gravity [5] . On the other hand, drifting populations are ubiquitous 
in complex plasmas. This occurs for ion-dust streaming instabilities [6] , elec- 
trostatic dust-cyclotron instabilities in collisional magnetoplasmas [7] and for 
current-driven ion-acoustic instabilities in complex plasmas [HI E] • 

In most of the above applications [H [10] , there is some streaming between 
electrons and/or ions relative to the dust, in the same direction of the mo- 
tion. In contrast, the present work consider a perpendicular drift of the ions, 
relative to the bulk dust motion. Ions and electrons are taken in local ther- 
modynamical equilibrium, which is equivalent to disregarding their inertia, 
while the dust fluid is assumed to be cold. Similar models were proposed 
for dust kinetic Alfven waves in dust-ion plasmas [TT| . but without the per- 
pendicular ion drift. As shown in the following, this ion drift is responsible 
for the generation of periodic magnetic field patterns. Notice that structures 
like braids, kinks and filaments at Jupiter's rings can be a consequence of 
an ambient magnetic field pattern and the associated force balance for the 
charged dust particles [T2] . 

This work is organized as follows. The necessary equations describing the 
electron, ion and dust fluids and the electromagnetic field are proposed in 
Section 2. Assuming a slab geometry, the system is then reduced to a nonlin- 
ear partial differential equation for the ^-component of the vector potential. 
This equation is studied in the linear and in the stationary limit. Linear 
solutions are shown to be unstable and to present a singular growth-rate for 
small wave-numbers. In addition, it is proved in the Appendix that the di- 
vergence of the growth-rate is related to a too fast characteristic time-scale, 
outside the range of validity of the model. On the other hand, the station- 
ary solutions are described by a nonlinear second-order ordinary differential 
equation, the Lienard equation, whose solutions shown periodic behavior. In 
Section III, the illustrative and relevant situation where electrons are fully 
depleted into the dust grains is discussed. Here the fundamental partial dif- 
ferential equation for the vector potential reduces to a Burger's like equation 
and there is no singularity in the dispersion relation. Section 4 is dedicated 
to the final remarks. 
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2 Basic equations 



Consider a plasma composed by electrons and ions, as well as by massive, 
negatively charged dust grains. As a first approximation, we take inertialess 
(in comparison to the dust particles) electrons and ions. For characteristic 
time scales much greater than the typical dust charge fluctuation time, the 
dust charge can be taken as a constant. Assuming a slab geometry, the basic 
equations are 

W + ^(n d v dx ) = 0, (1) 
= -Z d eE x - m d v dn v dx , (2) 
= e(Bl _ BoB J_^, (3) 

= - eEl -^k^L, (4) 
n e ox 

= ii e(ni - n i0 )v , (5) 

Z d n d = rn-n e . (6) 

In (PE]), the ion, electron and dust charges are e, — e and q d = —Z d e respec- 
tively, where Z d is the dust charge number and e the elementary charge. In 
addition, rii,n e and n d are the ion, electron and dust particle densities. The 
ions are assumed to have a perpendicular drift velocity —Voy, the negative 
sign chosen for later convenience. Drifting electrons could also be present, 
but here ions were preferred since in laboratory and astrophysical plasmas 
frequently the electron population is negligible (n e -C rii) due to electron 
charge depletion into dust grains [4j. Also, any current arising from dust 
grain flows is neglected, since it is smaller than the ion current by a factor 
(mj/m^) 1 / 2 at least, where and are the ion and dust masses, respec- 
tively. In this context, the magnetic field generation is due to the drifting 
ions only. Also, in the above equations B = B z (x,t)z is the magnetic field, 
E x the x component of the electric field, Vd x the x component of the dust 
velocity, v dn the dust-neutrals collision frequency and T i e the ion (electron) 
temperature, while the dust fluid is cold. Collisions between electrons and 
neutrals or between ions and neutrals are not taken into account. From 
Ampere's law (jSJ) it is implicit that rii is the total and n,o is the average ion 
number densities, so that dB z /dx — > as n, — > n i0 . A similar model [13] was 
recently proposed for a current-carrying electron-ion, non-dusty plasma. 
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In fl§]), quasineutrality is assumed, which is a good approximation for 
low-frequency phenomena. For the same reason, the displacement current is 
neglected, as well as the additional drifts arising from the generated mag- 
netic field. In this context, disregarding these Hall-type currents is the most 
delicate assumption, more throughly discussed in the Appendix. There it 
is shown that these perpendicular motions can be discarded in a sufficiently 
slow time-scale. The ion drag force due to momentum transfer between ions 
and charged dust grains is also not taken into account. This imply, in par- 
ticular, that the treatment is adequate for the co-rotating distance at the 
main (A — C) Saturn's rings, so that \v^q — v d a\ < vn, where and Vdo are 
the equilibrium ion and dust speeds and VTi is the thermal speed of ions. 
For outer rings, the plasma ions overtake the dust grains, while for smaller 
distances the ions are slower than the dust fluid. In both cases, signifi- 
cant momentum transfer takes place and more detailed models are needed 
PQ. However, the present work is basically concerned with the emergence 
of magnetic field patterns in highly collisional dusty plasmas, so that these 
effects are neglected. The model is still valid for the co-rotating distance at 
Saturn's rings (where x refers to the azimuthal direction) and for specially 
devised laboratory experiments. 

From now on we use the rescaled variables 

n e = n e /n i0 , ni = ni/n i0 , n d = Z d n d /n i0 , (7) 
E x = eE x /(K B Tik ) , B z = e\ Di cB z /(K B Ti) , (8) 
x = k Q x, i=u t, v dx = k v dx /u , (9) 



for 



v Z^jjlp? 



h = t , = -^r , (10) 

ADiC m d v dn c z 

where c is the speed of light, X^i = (eo K s^i/(^o e2 )) 1/ ' 2 the ion Debye length 
and ujpi = (nj e 2 /(mjeo) 1//2 the ion plasma frequency. As will be seen in the 
continuation, ko plays a major role as a maximal wave-number for linearly 
unstable modes. Also, it is implicit from (jUJ) that all wave-numbers are 
normalized to ko and that all frequencies and growth-rates are normalized to 
u . 

For simplicity of notation, the bars will be dropped Then, if the electro- 
static and vector potentials are denoted by (f> = <p(x,t) and A = A y (x,t)y, 
so that E x = —d(p/dx, B z = dA y /dx, then ([2JSD and (jSJ) gives 

n e = pexp(Tj0/T e ), (11) 
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Hi 

n d 

Vdx 



exp(— </ 
exp(— 4 

d(j)/dx . 



Ay) , 

A y ) -pexp(Tj0/T e ) 



(12) 
(13) 
(14) 



where p = n e0 /n iQ is the ratio between the equilibrium electron and ion 
densities. Notice that the y or z components of the model equations as well 
as Faraday's law are of no concern, since ([HE]) form a closed system. 
Eqs. (jSJ) and ffT2]) allows to express 



d 2 A 



while then (TlBl) gives 



f) 2 A ( T 

1 + 4^-pexp l-±L(A y + ln(l + 



d 2 A 



dx 2 



dx- 



V r)) 



(15) 



(16) 



Taking into account (114til6|) . all quantities can be written in terms of A y 
alone. The only remaining equation is (pQ), which leads to 



0_ 

dt 



d 2 A y 
dx 2 



pexp (-^(A y + \n(l + ^)) 



dx 2 



d 2 A v 1 
+ - 

dx 2 2 



'dAy 
dx 



+ ^exp(-^(^ + ln(l 



dx 2 

+ Ay + 

d 2 A 



(17) 



dx 2 



)) 



The resulting equation seems to be too complicated to be treated in full 
generality. In this Section, only the linear and the stationary limits will be 
considered. In Section III, the complete depletion case (p = 0) will also be 
shown to be amenable to some more detailed analysis. 

In the linear regime, (jTTI) reduces to 



d_ 

di 



d 2 A y pTi (d 2 A 



dx 2 



T, 



dx 2 



P) 



8P_ 

dx 2 



d 2 A y 
dx 2 



+ A, U 



(18) 



Assuming A y ~ exp[i(kx — ut)} then gives the dispersion relation 

p)(l-k 2 )k 2 



to 



k 2 -{pT % lT e ){\-k 2 



(19) 
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Taking into account that p < 1 because the dust grains are negatively 
charged, (IT3|) corresponds to instability (Im(uj) > 0) provided 

/ rp X-l/2 

k 2 _ < k 2 < 1 , where fc_ = 1 + -JL , (20) 



otherwise there is an evanescent mode. In addition, in terms of dimensional 
variables, (1201) reveals that the wave-number &o corresponds to marginal sta- 
bility, with k given by (TTUT) . Also, for p — >• 1, so that — >■ 0, the growth-rate 
becomes smaller. In other words, the present instability is intrinsically de- 
pendent on the dust, disappearing for electron-ion plasmas. 

For k —>■ k_, the growth-rate becomes singular. In several applications, 
one has k_ as a number not too far from unity. This corresponds to wave- 
lengths which are small but still much higher than the ion Debye wavelength. 
Therefore, a breakdown of the quasineutrality assumption is not a sufficient 
reason for ruling out the singularity as unimportant. Moreover, if the Pois- 
son equation is taken into account it can be shown that a (slightly modified) 
singularity is still present. Also notice that the singularity disappears when 
p^O. 

The clue to solve the problem is to notice that, when the growth-rate is 
too large, all quantities, and the magnetic field in particular, would change 
in a fast time-scale, in contradiction to the hypothesis of the model. This is 
a direct consequence of the neglect of the additional drifts in the y— direction 
due to the generated magnetic field (Hall effect). To prevent the emergence 
of a fast time-scale behavior, the wave-number has to be further bounded 
from bellow according to 

k > k m > k- , (21) 

where the full expression for k m is derived in the Appendix, see (141 p . 

Besides the small amplitude case, another relevant limit is the stationary 
one, so that dA y /dt = 0. In this case, the left-hand side of (IT7|) identically 
vanishes. The right-hand side, after integrated once assuming, for the sake 
of defmiteness, that dA y /dx and d 3 A y /dx 3 vanishes at x — 0, we derive, after 
some simple algebra, 

/ exp [-(TjT^Ay + ln(l + d 2 A y /dx 2 ))} \ _ 
{ 1 P 1 + d 2 A y /dx 2 J " U • 
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Defining B = dA y /dx, one of the equilibria is found when the first factor in 
is zero, 

d 2 B dB 



^ +B dF + s = °- < 23 » 

This is in the class of Lienard equations [2] and is the same equation found 
for the stationary magnetic field of a collisional current-carrying (non-dusty) 
plasma [13]. Observe that, by differentiation of (j!5p . the magnetic fields 
corresponding to fl23l) are associated to zero electric fields and dust velocities 
Vdx- In addition, from ( fTTT) the electron density would be constant, while 
from ( fT2lfT3"l) it is seen that the ion and dust populations appear modulated. 
The Lienard equation has an exact parametric solution [El HI], Du t here we 
content ourselves with a numeric solution, as in Figure (HI), showing a typical 
periodic behavior for B(0) = 0, B'(0) = 5. 

For p ^ 0, however, there is a second possibility for the stationary mag- 
netic field, where the second factor in (1221) vanishes. After differentiating, 
this can be shown to imply 



l + i^ + ¥} B T x +B) = - (24) 



Ti , d 2 B T 4 ,_dB 

1 + 

The further rescaling 

fi = (i + |r 1/2 5, £ = (i + |r i/2 x (25) 

converts ^2M into the same Lienard equation as before, 

d 2 B ^dB 



+ B— + B = 0. (26) 
dx z dx 

Despite the fact that both stationary states are described by the same (up 
to a rescaling) equation, there remains the conclusion that there are two 
possible nonlinear structures, provided p ^ 0. 



3 The complete depletion case (p = 0) 

It is well-known that complete depletion of the electronic population is not 
possible, since the minimum value of p 2 is the ratio between the electron 
and ion masses, assuming T e = Ti and point dust particles [E]- However, 
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assuming p = is a useful approximation. In such dust-ion 
plasmas, the wave dynamics is determined only by the inertia of the dust 
grains and the pressure of the ion fluid [16J. This is a relevant scenario, 
in particular, for laboratory experiments on dust-acoustic waves |3], certain 
stages of protostellar cloud collapse [10], the Saturn's F ring [16] and for 
Halley's comet [T7] . 

Postulating p = 0, together with integrating (fT7j) once, leads to 

dB _ d 2 B B 8B_ 

dt dx 2 dx 

This is a much more simpler diffusive Burger's like equation, in comparison 
with (TlTI) . It is formally the same as the equation for the magnetic field in a 
current-carrying electron-ion (non dusty) plasma [T3] . 

The results from Section 3 about linear waves and stationary states di- 
rectly translate to (1271) . For the linear waves, p = in (fT9l) gives uj = — k 2 ), 
which is an unstable mode for k 2 < 1. Here there is no singularity. Con- 
cerning the stationary states, they are described by a Lienard equation, as 
apparent either from ff27|) with time-independent B or from fl22|) with total 
electron depletion. Inverting the rescaling fflUj) . there follows the growth-rate 

Z d miLU pi , i/ x2 \2 r -2 



7=^ P -tf , (28) 
m d i/ dn V c / 

where K is the dimensional wave-number. 

The resemblance between (12~T1) and Burger's equation suggests the use 
of a Cole-Hopf transformation, which converts the Burger's equation into a 
diffusion equation. This is a convincing method to see that the dissipative 
term of Burger's equation prevents the formation of a shock. Therefore, it is 
illustrative to consider the Cole-Hopf transformation also for (12] 



B—j*. (29) 
q ox 



where q = q(x,t) is a function to be determined. Inserting (129]) into (J2j 
yields 

dtdq/dt-Pq/dJ^ n = Q) (3Q) 



9 

<9g 9 2 g 



9x 

or 

+ g(lng + r(t)), (31) 



9t dx 2 
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where T = T(t) is an arbitrary function of time only. However, this gauge 
function can be eliminated by 



q = e 9 q, (32) 

where g = g(t) satisfies the linear equation 

<7-<7 + r = 0, (33) 

so that q satisfies the same equation (l3Tj) as q, but with T = 0. Since q and q 
produces the same magnetic field through (12"9"|) . one can define T = without 
loss of generality to arrive at 

dq d 2 q , , , 

^ +q\nq. (34) 
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dt dx 

It would be a diffusion equation, except for the last (nonlinear) term. 

The stationary states where dip/dt = are associated to the solutions of 
a Hamiltonian system, 

d2q dV xr V 2 n 1 \ 

^ = -d^' V= ^ nq '2 ) - (35) 

From Figure [2] for the potential V, one conclude about the existence of stable 
periodic solutions around q — 1 and hence about the periodic structure of 
the stationary magnetic fields. These results applies directly also to the case 
of a current- carrying electron-ion plasma [13J. 



4 Conclusion 

In this work, a perpendicular (with respect to the main flow direction) drift- 
ing ion population was considered as a possible source of periodic magnetic 
structures in dusty plasmas. The basic equations reduce to a partial differ- 
ential equation of diffusive type for the z-component of the magnetic field. 
For small amplitude solutions, it was identified a linear instability, while the 
time-independent solutions are described by Lienard's equation. The peri- 
odic character of the solutions for the Lienard's equation can be associated to 
nonlinear magnetic field structures in space and laboratory dusty plasmas. 
The results are applicable, in particular, for highly collisional dusty plas- 
mas, where there is a noticeable momentum transfer due to dust-neutrals 
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collisions. Natural improvements are the inclusion of the ion and electron 
masses, as well as taking into account ion-dust or electron-dust momentum 
exchanges. Also a relevant open question is about the (nonlinear) stability 
of the predicted magnetic field structures. 
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Appendix 

Conditions of validity for neglecting the Hall effect 



A sensible condition for disregarding the Hall currents is given by 



\Vdx\ » \Vdy\ 



(36) 



where v dy is the y— component of the velocity of the dust fluid. Indeed, if 
( 136]) is satisfied, than the drift motion would not be an essential ingredient. 

Taking into account the y component of the force equation for the dust 
and E y = —dA y /dt as the y— component of the electric field, one get 



V d y 



Z d e dA y 
m d u dn dt 



(37) 



where dimensional variables are used. It is evident from (13611371) that the 
magnetic field should not change in a too fast time-scale. While the condition 
is clearly achieved in equilibrium, this is not what happens in the linear 
regime when approaching the singularity k — > k_, with k- given in (I22|) . 
More exactly, using (|14p and fl37|) and from now on the rescaled variables, 
f l36|) is rephrased as 



OA. 



dt 



« 



m d cv dn 



ZdlJliVTiUJ' 



pi 



dx 



(38) 



In ([3]), vt% = (KBTi/rrii) 1 ^ 2 is the ion thermal velocity and the ion plasma 
frequency is uj p i = (njoe 2 /(mjeo)) 1 ^ 2 . Except if the electrostatic potential 
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significantly changes in a short spatial distance, (158)) is usually fairly satisfied 
both in laboratory and astrophysical plasmas. 

More detailed conclusions can be found in the linear regime, which is also 
the most problematic. Indeed, diffusion and nonlinear effects are expected to 
prevents the development of the singularity predicted for k — ► fc_ . For small 
amplitude waves, (ITB1) gives 



Taking into account (1391) and after Fourier transforming (1381) . there follows 

l m {u) « mdCUdn k \l - k 2 \ , (40) 



Z d miV Ti uj, 



pi 



where uo is given by the dispersion relation ( fT9l) . Now some algebra shows 
that (PR)]) is equivalent to k » k m , where 



P + 



1/2 

'1 n "\2 _i_ A Tjp l m d cu d „ \2/i I Tip\ 

v 1 W I - * T e \Z d m lVTl u v J \ L ^~ T e > 



2m d cv dn (l + ^)/(Z d miV Ti u pi j 



(41) 



In particular, for the fully depleted electrons case, the smallest allowable 
wave-number reduces to the much simpler form 



Z d miV T iU) p 



* pi 



m d cv dn 

while, in the absence of dust, it reduces to 



p~0, (42) 



/ T \- 1/2 

km= \ 1 + ¥J ' p ^ 1 ' (43) 

For instance [7j, in comet tails with n e o — n i0 = 10 cm" 3 and a small fraction 
n d o/n i0 ~ 10~ 6 , one can apply (1131) with T e /Tj ~ 5, yielding k m = 0.41. Also 
notice that in the pure electron-ion case, k m as given by (1431) is the same as 
k_ in (120|). 

As apparent from either (j4"Tl) or (pj2]), without discussing the role of the 
remaining parameters, the range of validity of the model increases for a 
larger dust-neutrals collision frequency. This is natural, since perpendicular 
drifts tends to be inhibited in a strongly collisional medium. As an example, 
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consider the plasma parameters in the vicinity of Saturn's main A and B 
rings, as recently measured by the Cassini spacecraft orbiter and shown in 
Table 1 of pQ. Such planetary rings are basically a dusty plasma with low 
degree of ionization. Indeed, the Cassini measurements reports on an oxygen 
atmosphere and an ionospheric plasma near Saturn's main rings. There the 
ratio between the ion and neutral molecular oxygen 2 densities is as 
low as 10~ 3 — 10 -4 , so that ion-neutral momentum transfer is expected to 
be a relevant process. For [T] the plasma parameters p ~ 1/10, Tj ~ T e ~ 

1 - 10 eV, m d = 4 x 10- 12 g, v dn ~ lO" 11 - lO" 10 ^ 1 , Z d ~ 1 - 10 2 , v Ti (Ot) ^ 

2 x 10 5 - 7 x 10 5 cms' 1 , u pi (0£) ~ 2.5 x 10 2 - 1.2 x 10 3 s -1 , one get from 
(l4"Tj) the range k m ~ 0.3 — 0.7. Better results could be found for outer 
rings like Saturn's F ring, where p has much smaller values |4]. For these 
distances, however, the plasma ions significantly overtake the dust grains, 
so that momentum exchanges between ions and charged dust grains should 
be considered pp. In all cases, the singularity predicted by the linear theory 
for k = k_ is unphysical, since /c_ as given by (120]) is smaller than k m . 
Considering the large degree of depletion, if one uses fT2"8l with K = and 
Vq/c = 10~ 2 one roughly estimate the maximum growth-rate as 0.8 — 1.9 x 

ioV 1 . 

A practical way to reach smaller wave-numbers is to enlarge the dust- 
neutrals collision frequency, which can be calculated [7] from 

Am N Na 2 v TN 

v dn ^ , (44) 

m d 

where is the neutrals mass, iV the neutrals number density, vtn the 
neutrals thermal speed and where the dust grains where assumed to be all 
spheres of radius a. Therefore, for space plasmas or for laboratory experi- 
ments, for smaller fc m 's one should have, in particular, the largest possible 
neutrals mass density combined to the largest neutrals thermal velocity. In 
addition, it follows from fj4T|) that larger wave-lengths are achieved for larger 
values of T e /Ti. 

As a final example, one can choose parameters appropriated to the dust- 
acoustic instability. The dust-acoustic experiments are promising, since they 
consider a small ionization fraction (~ 10 -5 ), so that the ion drag force 
is smaller than the neutral friction even if v dn is smaller than the ion-dust 
collision frequency 0. Following [9], take p = 1/2, Tj = 300K, T e = 23000X, 
mi = m N = 4.7 x 10~ 23 g, m d = 10 _9 g, Z d = 4 x 10 4 , v Ti = v TN = 3 x 
10 4 cms _1 , u P i = 10 3 s _1 , iV = 3xl0 15 cm~ 3 , a = 5 x 10~ 4 cm. Combining fT4TT) 
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and (j44j) . there follows k m = 0.08. This small value comes basically from the 
large ratio T e /Tj. Finally, the maximal growth-rate can be estimated using 
(j!9p . since now the depletion is not so high. Considering (fT9l) with k = k m , 
and for vq/c = 10~ 2 , there follows the value 8.1 x 10 3 s _1 for the growth-rate. 
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